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Abstract 

The Lieb-Robinson theorem states that locality is approximately preserved in 
the dynamics of quantum lattice systems. Whenever one has finite-dimensional 
constituents, observables evolving in time under a local Hamiltonian will essen- 
tially grow linearly in their support, up to exponentially suppressed corrections. 
In this work, we formulate Lieb-Robinson bounds for general harmonic systems 
on general lattices, for which the constituents are infinite-dimensional, as systems 
representing discrete versions of free fields or the harmonic approximation to the 
Bose-Hubbard model. We consider both local interactions as well as infinite- 
ranged interactions, showing how corrections to locality are inherited from the 
locality of the Hamiltonian: Local interactions result in stronger than exponen- 
tially suppressed corrections, while non-local algebraic interactions result in alge- 
braic suppression. We derive bounds for canonical operators, Weyl operators and 
outline generalization to arbitrary operators. As an example, we discuss the Klein- 
Gordon field, and see how the approximate locality in the lattice model becomes 
the exact causality in the field limit. We discuss the applicability of these results 
to quenched lattice systems far from equilibrium, and the dynamics of quantum 
phase transitions. 



1 Introduction 

Locality in relativistic theories ensures that space-like separated observables commute: 
One simply cannot communicate faster than light. In non-relativistic lattice models, 
in contrast, there is no a priori reason for the support of time-evolved operators to 
stay confined within a light cone. It is one of the classic results of mathematical 
physics, dating back to Lieb and Robinson JH], that even in a non-relativistic quan- 
tum spin model on a lattice, locality is preserved in an approximate sense: There is 
always a well-defined velocity of information propagation and hence a causal, "sound" 
or "light" cone. Locality is then respected under quantum spin dynamics with finite- 
ranged interactions, in that the support of any local observable evolved for some time 
will remain local to a region of size linear in this time, up to a correction that is at 
least exponentially suppressed. Except from exponentially decaying tails, hence, one 
encounters a situation very much like in the relativistic setting. The Lieb-Robinson 
theorem has hence put the observation in many quantum lattice models on a rigorous 
footing that there exists a well-defined finite speed of propagation, often referred to as 
group velocity. 
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Over the years, Lieb-Robinson bounds have been extended and generalized to 
higher-dimensional spin systems on lattices, and the bounds have been significantly 
improved in several ways llJ|2l|3j|4ls|aLZ|,|8|]. Also, important new applications of 
the Lieb-Robinson theorem have been found: Notably, the only known proof for the 
clustering of correlations in gapped lattice models is based on the Lieb-Robinson theo- 
rem hence rigorously confirming the "folk theorem" in condensed matter physics 
that correlation functions decay exponentially in gapped models. "Area laws" for the 
scaling of entropy in grou nd states of one-dimensional spin systems can be proven 
based on this result lf9L 1 1 OL 11]. Finally, in the context of quantum information theory, 



it provides a bound to the velocity one can transmit quantum information through a 



chain of systems giving rise to a quantum channel 111 111 , a topic that has received quite 
some attention in the quantum information literature. 

In the simplest form of the Lieb-Robinson theorem, one considers a spin system 
H = ©ieL %> % = ^ > on a lattice with vertices L, and a local (finite-ranged) 
Hamiltonian H = J^xcl ^(^0- The time evolution of an observable A on some 
subset A C L of the lattice under this Hamiltonianj]] 

A(t) = e i6t Ae~ iAt , (2) 

then forms a group of automorphisms. At time zero, obviously any observable A being 
supported on A will commute with any observable B that is supported on a disjoint 
set B C L. The Lieb-Robinson bound now gives a bound to this commutator if A is 
evolved in time under this local Hamiltonian. It says that there exists a constant C > 
and a "speed of light" v > such that 

\\[A(t),B]\\ < C\\A\\ \\B\\c-^ d ' lst{A ^- vW \ (3) 

where dist(A, B) = m\n a£ A.b£ b dist(a, b) is the minimal distance between the two 
regions and || • || denotes the operator norm. In other words, outside the causal cone 
v\t\ < dist(A, B), one encounters merely an exponential tail, and the supports of 
A(t) and B stay almost disjoint. Eq. © governs the maximum speed at which a local 
excitation can travel through the lattice and the maximum speed at which correlations 
can build up over time. 

The physical setting considered here can equally be viewed as the study of the sit- 
uation of quickly quenching from, say, a system that is in the ground state of some 
local Hamiltonian to a new local Hamiltonian 



setting has also been linked to the entanglement generation and scaling in quenched 
systems Hoi [Til IT2I [Lli l23ll . Studies of non-equilibrium dynamics of quantum lattice 
systems of this type are entering a renaissance recently, not the least due to experi- 
mental studies becoming more and more available. With atoms in optical lattices, for 
example, one can suddenly alter the system parameters and thus observe the dynamics 



1 This automorphism group of time evolution is in the context of Lieb Robinson bounds typically denoted 

as 

r t H (A) = c iHt Ac- iHt . (1) 

Since in this work, the Hamiltonian H is always well-specified at the beginning of each section, we make 
use of the above notation in the Heisenberg picture for simplicity of notation. 
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of a quantum phase transition 12411 . Hence, it seems only natural to apply the machin- 
ery of Lieb-Robinson bounds to such settings. However, despite the generality of the 
above mentioned results on Lieb-Robinson bounds, they are — with the exception of 
the recent Refs. |0] — constrained in the sense that they only apply to spin systems, so 
finite-dimensional constituents, quite unlike the situation encountered in many settings 
of non-equilibrium dynamics. 

In this work, we derive Lieb-Robinson bounds to harmonic lattice systems on gen- 
eral graphs. Such models correspond to discrete versions of free fields, lattice vi- 
brations, or the superfluid limit of the Bose-Hubbard model. As harmonic models, 
applicable merely to a class of systems, the resulting bounds can indeed be made very 
tight, e.g., for local Hamiltonians, we find stronger than exponential suppression, while 
for algebraically decaying interactions the corrections to locality are also algebraically 
suppressed. Within the considered class of models the very tight connection between 
the approximate locality of operators and the locality of the Hamiltonian is thus re- 
vealed. 

These systems serve as instructive theoretical laboratories for more elaborate inter- 
acting theories for infinite-dimensional quantum systems (and Lieb-Robinson bounds 
have fundamental im plic ations, e.g., in the context of using harmonic systems as quan- 
tum channels 1 12 , 13lll4ll ). about which little is known when it comes to Lieb-Robinson 
bounds (see, however, Refs. 10,01 for recent prog ress on bounded anharmonicities). In 
this way, we continue the program of Refs. l25l 12^1 . building upon earlier work pri- 
marily on clusteringof correlations and "area theorems " in harmonic lattice systems 

iSlilSSEllii. 

This chapter is organized as follows. We first define the models under considera- 
tion, and explain what we mean by having a general lattice. We then present bounds 
on the time evolved canonical coordinates and make the causal cone explicit. For the 
important class of Weyl operators on the lattice, we explicitly find bounds on their op- 
erator norms and discuss generalizations to arbitrary operators. As an example, we 
discuss the case of the discrete version of the Klein-Gordon field, and show how the 
approximate locality in the lattice model becomes the exact locality in the continuum 
limit. The proofs of the main results are then presented in a separate section in great 
detail. As such, the findings in this work complement the findings of Ref. jg], which 
considers in its harmonic part nearest-neighbor interactions of translationally invariant 
models on cubic lattices. We finally discuss the implications to quenched many-body 
systems far from equilibrium. 



2 Considered models and main results 

We consider harmonic systems on general lattices. Such lattices are described by an 
undirected graph G = (L, E), with vertices L and edge set E. The vertices L cor- 
respond to the physical degrees of freedom, here bosonic modes with Hilbert space 
Hi = £ 2 (R), i € L. Edges reflect a neighborhood relation on the lattice. On L we 
use the graph-theoretical distance dist(i,j) between i, j € L, i.e., the shortest path 
connecting i and j. On such a general set of lattice sites L we consider harmonic 
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Hamiltonians of the form 



H = g ftiXijXj ■ /',/', ,/',! , (4) 

where Xjj = X,-^ e R, Pi.j = Pj t i € R and the aij, p, are canonical coordinates 
obeying the usual commutation relations (we set H = 1) = [pi,Pj] = 0, 

[xupj] = iSij. Identifying 

a _ + Pj.j R _ Xjj — Pj j - _ Xj + ipi 

the above Hamiltonian is equivalent to a Hamiltonian quadratic in the annihilation and 
creation operators of the bosonic modes ([Si, bj] = 0, [hi, &j] = Sij) of the form 

6 = \ £ ( h l A iJi + kAj] + biBiJbj + blB itj b]) . (6) 

All the following results may thus also be obtained for the above Hamiltonian and 
bosonic operators by using the identification in Eq. ©. We suppose that L is countable 
such that we may identify the couplings [Ci,j)ije.L> C = X, P, with matrices, thereby 
defining the operator norm ||C|| and multiplications in the usual matrix sense. We 
denote the time evolution of operators 6 in the Heisenberg picture as 

S(t) =e lkt oe-' lk \ (7) 
and by 1 1 ■ 1 1 the operator norm throughout. 

2.1 Local couplings 

In this subsection, we will derive Lieb-Robinson bounds for harmonic systems with 
arbitrary local interactions on the graph. We will see that outside a causal cone we 
obtain a stronger than exponentially suppressed influence of time evolved canonical 
coordinates in the Heisenberg picture. In the above notation, local means that 

Xi j = Pi j = for dist(i, j) > R. (8) 

For notational clarity we write 



dij := dist(i,j)/R, t := max{y/\\PX\\, y/\\XP\\}\t\ (9) 

in the following, and denote by \x] = min{z € % \ x < z} the ceiling function. 
For a Hamiltonian as in Eq. © with local couplings as above we prove the following 
theorems. 

Theorem 1 (Lieb-Robinson bounds for local couplings) Writing b^j = \d h j/2\ and 
a,j = max{0, |~(djj — l)/2]}, one has 

[ x i(t),Xj]\\ , uvu II lPt(t),Pj}\\ < 



P\\ " L JJ "' ||X|| (2a,, +1)! ' (1Q) 

\\Mt),p-M,\\lMt)^M< T 



(26,,,) 



4 



We note that (d€N and we use l/d\ < (e/ d) d d~ 1/2 ) 

T d coshM e r+d ( 1+log ( r )~ log ( d ^ 

T -^<- 73 • 

i.e., for sufficiently large dist(i,j), one finds a faster-than-exponential decay. In the 
subsequent formulation we make this more explicit by defining a "light cone", er < 
dij (i.e., c\t\ < dist(i,j)), with a "speed of light" given by 

c := ci?max{||XP|| 1/2 , ||PX|| 1/2 }. (12) 

Then commutators of "space-like separated" operators are strongly suppressed. This c 
is an upper bound to the speed with which a local excitation would travel through the 



lattice in a non-equilibrium situation. This kind of argument is used, e.g., in Ref. II15H . 
where a central limit type argument was used to show exact relaxation in a quenched 
system - the intuition being that inside the cone excitations randomize the system while 
the influence of excitations outside the cone is negligible, which is essentially a Lieb- 
Robinson-type argument. In a very similar fashion, one can argue that the speed at 



which correlations build up in time is governed by these bounds (see also Ref. 11 111 ). 
Again, the above speed is an upper bound to which one can correlate separate regions 
starting from an uncorrected state under quenched, non-equilibrium dynamics. 

As mentioned before, these bounds have immediate implication to the evaluation of 
capacities of harmonic chains, when being used as convenient quantum channels. Such 
an idea of transporting quantum information through interacting quantum systems is 
an appealing one, as no exact local control is required, and the transport of quantum 
information is merely due to the free evolution of an excitation under local dynamics^ 



Following the argument of Ref. [11], one could use such a harmonic chain as a quantum 
channel through which one sends classical information, encoded in the application of 
a local unitary at some site, and letting the system freely evolve in time. The decoding 
corresponds to a readout at a distant site. Then indeed, outside the cone defined by 
the "speed of light", the classical information capacity C of this quantum channel is 
exponentially small. This means that the classical information capacity of harmonic 
chains used as quantum channels is - for a fixed time - exponentially small in the 
distance between sender and observer. The causal cone is made even more explicit in 
the subsequent formulation. 

Theorem 2 (Alternate version making the causal cone explicit) Leter < dt.j. Then 

IU Xi(t),Xj]\\ , ^TViP ll[Pi(*)>#/]ll J II > and II II 



||P|| ' V " JJ "' ||* 

are all bounded from above by 



di,j log(eT/di,j) 



(i - (er/^,i) : 



(13) 



(14) 



^For harmonic instances, see, e.e., Refs. (JHHGl, 

but there is a vast literature also for spin systems 
and other finite-dimensional quantum systems, to mention all of which would be beyond the scope of this 
chapter. 
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Often, X and P will commute, rendering the max in Eqs. ( 19112b irrelevant. In 
many physical situations one even has Pij = Si_j, i.e., r = \J \ \ X j 1 1 1 | , in which case 
the above bounds may be improved, in particular, the "speed of light" improves to 
c = ci?||Xj| 1 / 2 /2. For clarity, we explicitly state the new bounds in this case, in both 
ways. 

Theorem 3 (Lieb-Robinson bounds for local couplings and P = 1) Letaij = max{0, \ckj- 
1]}. Then 

T a r*,/l +i cosh (t) 



_ii b , (t) , Pjlll <__U, (15) 

r 2r*,jl cosh(r) 



ll[^(<),^]IMI[^W,%]ll < 



(2^-1)! 



Theorem 4 (Alternate version making the causal cone explicit for P = 1) For cr < 



2<ij ] j one /ifli 



and for ct < (2ajj + 1) 



Jdij log(eT/(2d i , j )) 



/d~ (l - (er/(2d iji )) 2 ) 



(16) 



I „2a i ,3log(eT/(2a i , j +l)) 

"^(i),^-]!! < t ^ =v, (17) 



' y/aij (l - (cr/(2a t . J + 1)) 

where now a,ij = max{0, \di.j — 1]}. 

2.2 Application: Non-relativistic quantum mechanics yields causal- 
ity in the field limit 

This section forms an application of the previous considerations. We will see how the 
exact light cone of the free field is recovered from the approximate light cone in the 
Lieb-Robinson theorem in the continuum limit of the lattice version of the field theory. 
It is very instructive indeed to see how the tails in the superexponentially suppressed 
region outside the light cone becomes more and more suppressed in this limit. The role 
of the Lieb-Robinson velocity is hence taken over by the speed of light in the relativistic 
model. 

We start from the Klein-Gordon Hamtttonian on V = [0, 1] xD in units fl = c = 1, 



H = - [ dx 

2Jv 



(18) 
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where the field operators fulfill the usual commutation relations 

[<p(x),7t(y)]=i6(x-y), [<p(x), <p{y)] = [w(x), n(y)} = 0. (19) 
Discretizing according to x = i/N, i G {j G 1N D \ id = 1, 2, iV} =: L, 



I/AT 



(20) 



where denotes a unit vector in direction d, we find (equipping L with periodic 
boundary conditions), 



H 



— Y 

2N D 

i£L 
( 



1 



2N D 
='■ Hn. 



^(i/N) + J2 

d=l 



'<p(i/N + n d /N)-tp(i/N) s 2 



1/N 



m 2 ip 2 (i/N) 



\ 



2DN 2 )<p 2 (i/N)]-N 2 Hi/N)0(j/N) 



i.jeL 

dist(i,j) — l 



J 



(21) 



Then N — > oo is the valid continuum limit for a fixed V — [0, 1] xD . Now, 

&i := N- D ' 2 (p{i/N), pi := N- D / 2 n(i/N), (22) 

define harmonic position and momentum operators satisfying the canonical commuta- 
tion relations, in terms of which we find 



H n — 



ijeL 



where Pij — Sij and 



Xi 



(m 2 + 2DN 2 )S i j-N 2 5 dlst{iJhl . 



(23) 



(24) 



We are interested in the discretized version of the commutator [<p(x, t), <p(0, 0)], which 
is given by 



N [xi(t),x ], 



(25) 



and set out to apply Theorem |4] We have R = 1 and assume w.l.o.g. that < id < 
N/2, i.e., = dist(i,0) = J2d=i^ d = ^J2d=i Xd — N\x\, with |x| being the 
euclidean norm. Furthermore, as we assumed translational invariance the eigenvalues 
Afc of X are given by 



D 



A fe = m 2 + 2DN 2 - 2N 2 ^ cos (2irk d /N) , 



(26) 



<2=1 
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i.e., \\X\\ < m 2 + 4DN 2 (for even N we have equality). 
Now fix It I and \x\ such that 



eVB\t\ < \x\. (27) 



We then take the limit N — > oo such that J2d=i Xd ~ di,o/N — const, is fulfilled for 
all N (e.g., x = (1/4, 0, . . . , 0) fixes i = (N/4, 0, . . . , 0) and N/4 € N). Now let 
N G N be such that 

e\t\VD I m 2 



which yields 

er < | g | vW < 2rf . < 2di0 (29) 

for all > Nq. This enables us to apply Theorem|4]to find 

2Ar|x|log( Z ) 

Pi(*),£o]|| < (30) 

v^ll^lll^l (i - z ) 

for all N > N , i.e., 

lim N D \\[xi(t),x ]\\ =0 (31) 

independent of m. Eq. (1311 shows that the approximate light cone of the Lieb-Robinson 
bound becomes an exact light cone in the continuum limit. The exponentially sup- 
pressed tails vanish, and approximate locality is replaced by an exact locality. It is 
interesting to see how this concept emerges from the bounds to the speed of informa- 
tion propagation in the sense of Lieb-Robinson. 

The bound in Theorem |4] is not quite strong enough to recover the exact prefactor 
of the light cone \t\ < \x\. This is mainly due to the fact that we allowed for gen- 
eral lattices in the Lieb-Robinson bound. Demanding translational invariance would 
allow for slightly stronger bounds. In Fig. Q] we depict exact numerical results for this 
geometrical setting in D = 1. 

2.3 Non-local couplings 

The previous section allowed for arbitrary local interactions. In this section we will 
turn to strongly decaying non-local couplings of the form 



(dist(i,j) + 1) 



\X id \ , \Pij | < 717TIFl \ ; ^ . (32) 



We define the spatial dimension of L in the usual sense: For all spheres S r (i) C L with 
radius r G IN centered at i £ L, 

S r (i) = {ke L\dist(k,i) = r} , (33) 
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N\\[xi(t),xo]\ 




f 



250 



N 
500 




*IIN*-IBo).*o]|| 



N 



Figure 1: The light cone in the field limit of the discrete Klein-Gordon field: De- 
picted is N\\ [£i(t), xq] I as a funtion of t and N. This is the discrete version of 
|| [<p(x, t), ip(Q, 0)]||, x = i/N, in one spatial dimension. For \t\ < \x\ we find for 
finite system size N the Lieb-Robinson exponential decay in \x\ and for N — > oo 
the commutator N\\ [xi(t), xq]\\ goes to zero for space-like separations, thus recovering 
exact causality. 



there exists a smallest D > such that for all < r £ N 



sup |<SV(*)| < cd?" 

i£L 



D-l 



(34) 



for some constant cd > 0. This number D is taken as the dimension of the lattice. 
We find that the decay of interactions is inherited by the decay of the operator norm 
of the commutator of canonical coordinates. The same power in the exponent as in 
the interaction again appears in the Lieb-Robinson bounds. Note the (not accidental) 
similarity with the inheritance of the decay of correlation functions dependent on the 
decay of interactions in Ref. 025I1 . 



Theorem 5 (Bounds for non-local couplings) Let > D. Then 

sinh(r) 



\\Mt)^i]\\MMM\\ < 

\\[xi{t),pj\\\,\\]pi(t),Xj]\\ <5 



a (l + dist(i,j)yi' 
cosh(r) 



(35) 



h3 a a (l + dist(t,j))^ 
where t = aoCo\t\, clq = C£>2 ?7+1 C(1 — D + 77), and £ is the Riemann zeta function. 

2.4 Weyl operators 

A class of operators that play a central role in harmonic systems are the Weyl operators. 
Denoting the support of a Weyl operator as S C L it may be written as 



where£= (xi, zisbPi, — ,P|s|) € IR. 



(36) 



9 



Via the Fourier- Weyl relation general bounded operators may be expressed in terms of 
these operators, see below. 

We define the distance of two sets A, B C L as 



dist(A,B)= min dist(i,j), 

and the surface area of a set A C L as \dA\, where 

dA = {i e A\3j e L\A : dist(i,j) = 1} 



(37) 



(38) 



defines the set of lattice sites on the surface of A. The following theorem establishes a 
connection between commutators of Weyl operators and previously derived bounds on 
the canonical coordinates. Note in the subsequent theorem the dependence on the right 
hand side of the operator norms 1 1£| | , 1 1£' 1 1 of £ and whereas on the left hand side of 
Eq. (1401 . we have the operator norm for commutators of Weyl operators. 



Theorem 6 (Lieb-Robinson bounds for Weyl operators) Let 



be Weyl operators as defined above. Then 



(39) 



<u\\\\a E (ii[^(*).^]ii + ii[^(*)»fi]ii 

+ llbi(*)»^]ll + llbi(*),Pj]ll) 

< CD u\\u'\\^{m\,m'\} 

oo 

x f(d)d°- 1 (l + c D (d-dt S t(E,E')) D ) . 

d—dist^^') 

where f : N — > 1R is a function such that 

\\[xi(t),*j]\\ + \\[xi{t),Pj]\\ + \\\Pi(t),Xj]\\ + WlMt),PM < f(dist(i,j)). 
Employing, e.g., Theorem^ we have for er < dist(S, =: ds.H' that 



(40) 



(41) 



Wt(t),W e \ <Cmm{\dZ\,\dE'\}g 
where 

C = R^cnUWU'W 
and the function g : (0, 1) — > R, 



3 d H _ H , lo S ( e T/d BiB ,) d D-3/2 



\X\ 



is increasing in z with lim 2 



z d ' R (d + 1) D ~ X {1 + c D (d+ l) D ) > 0, 



(42) 



(43) 



(44) 



d=0 
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Note that we have in Eq. (HTt expressed this statement in terms of a function / that 
grasps the decay of operator norms of commutators of canonical coordinates. When- 
ever one can identify such a function, e.g., through theorems 1-5, a result on Weyl 
operators can be deduced. Needless to say, in the same way we have applied Theorem 
12 we could have made use of Theorem [5] Essentially, the decay in operator norms 
of canonical coordinates is inherited by the expression for Weyl operators. Due to 
the sum in Eq. d40t . however, it can happen that no decay follows for Weyl operators 
if the (i) dimension of the lattice is too high or (ii) the decay of / is too slow. For 
finite-dimensional lattices, however, sufficiently fast algebraically decaying (i.e., suf- 
ficiently large if) interactions yield an algebraic Lieb-Robinson-type statement for the 
commutator of two Weyl operators. 

Note also that only the surface areas of the two sets S, S' enter the bound, but not 
the cardinality of the support. This allows for infinite regions (for D = 1 both may 
be supported on infinite intervals, for D > 1 only one of S, 5' needs to have a finite 
surface area) separated by dist^E, S'). 



2.5 More general operators 

A general bounded operator 6 supported onScL may be expressed as 



1 



where 



(2tt)W 



Xo(0 



]R 2 I H I 



oW f 



(45) 
(46) 



is the characteristic function of o. This allows to deduce bounds for general bounded 
operators using the bounds on Weyl operators stated above: 

f d£ \ X o(-0x6'(-e)\ \Wz(t),W?] . (47) 



m,o'}\< 



R 2 I S I </]R 2 l s/ 



(2tt)I s I+I h 'I 

Bounds for more general, possibly unbounded, operators that are finite sums of 
finite products of canonical operators (or bosonic creation and annihilation operators) 
may be obtained by repeatedly employing operator identities such as 

(48) 



Ab,c 


= A 


B,C 


+ 


A,c 


B, 


A,BC 




A,B 


C + B 


A,C 



e.g., bosonic density-density commutators may be written as 

[Ai(t),ftj] = fet(t) [8i(*),ftj| + [&!(*), k(t) 



Hit) (\k(t)M 



bi(t),bj 
h + b) 



(49) 



b}(t),bj )bi(t), 



which yields 

mm,^ < (b\(t)b 3 



k(t),b] 



b\{t)b\) k( t ),b 



b]k(t) 



W),b 3 



(50) 
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where bounds on the commutators may then be obtained by identifying bosonic oper- 
ators by canonical operators through Eq. © and employing the above derived bounds. 

3 Proofs 

In this section, we will present in detail the proofs of the previous statements. 

3.1 Preliminaries 

We write the Hamiltonian in Eq. © as 

1 2|i| 

H= 2 X! ^n,H nunj f nj , (51) 

ni,rirj = l 

where we have arranged lattice sites such that H ni , n = Xi,j (= Pi,j ) for 1 < n,, rij < 
\L\ (L+l < n^n.j < 2| J L|)andf„ i = x { (= ft) fori < n- < \L\ {L+l < m < 2\L\). 
Now consider the time evolution of the operator 

f ni (t)-=e i6t f ni e- i6t . (52) 

By solving Heisenberg's equation of motion or, alternatively, by employing the Baker- 
Hausdorff formula, one finds 

2\L\ 

r nt (t) = ^2 ( e ~ aHt ) ni . nj where Un ^ n i = 1 [fni,f ni ] ■ (53) 

nj — 1 

This yields for the commutator 

2\L\ 

i [f ni (t),r nj ] = i £ (e ""•')„ [r„ fc ,f„J = (e-° Ht *) nunj . (54) 

Now, separating the terms with an even power in n from the terms with an odd power, 
we get 




~ (_l)n f 2n / (p X) n Q \ 

h ( 2 ") ! V o {xpf ) ■ 
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Hence, 

1 = E (2n+l)! = : 1 ' C ^ 

n=0 ^ ' ' 

00 / 1 \nj-2n 



(56) 



n=0 



These expressions will form the starting point of the subsequent considerations. 

3.2 Local couplings 

We will need the following lemma. It states that finite powers of local coupling matrices 
defined on graphs remain local couplings, albeit with a larger range. 

Lemma 1 (Products of local couplings) Let A = (Aij)i^,=L be such that Aij = 
for dist(i, j) > R. Then for n G N 

(A n )i t j = Ofor all i,j G L with dist(i,j) > nR. (57) 

Proof. For 11 = 1 the statement is obviously true. Now let (A n )ij = for dist(i,j) > 
nR. Then 

Now let k G L. If dist(k,j) > R then this k does not contribute to the sum as 
Akj = 0. Now let dist(i,j) > (n + l)R. Then we have that also if dist(k,j) < R it 
does not contribute to the sum as then dist(i, k) > nR and therefore (A n )ij = 0: 

(n + 1)R < dist(i,j) < dist(i, k) + dist(k, j) < dist(i, k) + R, (59) 

i.e, dist(i, k) > nR. □ 



Thus, if we have Xij = Pij = for dist(i,j) > R, we may write (see Eqs. (l56ll), 

|y-l 2 ™+ 1 

(60) 



00 U|2n+1 

\ c m\ < E ^n-TTy.\ iiPxrp) 



>-,j 



where ctfj = max{0, \{<U,j ~ 1 )/2~| } and we recall that dij = dist(i, j)/ R. As one 
has for any matrix that \ < \\M\\, we find 

llpll ~ r 2n 
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where we recall that r = max{ i/||PZj| , y^jjXPj[}|*|. Similarly (b id := \d itj /2\) 



|c^(*)|,|c^(*)|< E 



(62) 



(2n)! 



In the case of = djj these bounds read 



I 00 T 2n 

l^jWl ^ /TmFTT E To^m 



vW«=r«^ + i/ a (2n)! 

00 2n 

|CS(*)|<>/Pil E T^v (63) 



00 2n 

ra*)i.ra*)i< e 7b 



(2»)! 



Now, for c > 0, 



r 2c r 2n T 2c 



g (2^)! = E (o^)! ^ n ?„W = M C ° Sh(T) ' ^ 
also, if er < 2c, we find 

f Z!l < V (2n)-/2 < J_ f (*_LY n = 

^ c (2n)! ~ ^\2n) V 7 " ^Tc ^V2cV ^ (1 - (er/2c) 2 ) ' 

(65) 

where we have used that n! > (n/e) n n 1 ^ 2 for n > 1. 

3.3 Non-local couplings 

Let \Mij\ < [1 + dist(i, j)]~ v - For such couplings we have (dij := dist(i,j)) 



'M 2 ) 



k 

/ 2 \" 

< 



2 ^ n ( 1 + max{d i)fc , } ^ ^ 



1 + <W ^ V(! + rfi.fc)(l + rffe.j) 



where we have used the triangle inequality and (a + b) < 2 max{a, b}. Now, 

(1 + max{d i , 1e ,d kJ }) r i < (1 + dj.fc)" + (1 + (67) 
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i.e., an upper bound for the above sum over k is given by 
1 1 



E 



(l + d t . k )v (l + dkj) 71 



< 



= *E 



(l + r) 7 ' 



where 



|{fcGX|d M =r}| = |5 r (i)|, 



which we may bound using the definition of the dimension of the graph to find 

2^+i °° i 



A/ 



which converges if r\ > D, in which case we have 



M 



< 



no 



no 



where £ is the Riemann zeta function. By induction we then find 



< 



-o"o 



(1 + 



(68) 



(69) 



(70) 



(71) 



(72) 



for n > 1, implying (recalling that r = coao|t|) 
1 



\cfm\ < 



)n E 



r 2n+l 



sinh(r) 



o (l + dist(i,j)) r > ^ (2n+ 1)! o (l + dist(i, j)) 7 * ' 



\C%(t)\<6 itj 



and similarly 



1 



a (l + dist(i,j))v ^{2n)\ 



\v E 



00 2n 



= *ij + 



cosh(r) 



(73) 



a (l + dist(i,j))^'' 



\cm\ < 



sinh(r) 



^ y y| " aoil + distiij))*' 
cosh(r) 



(74) 



\C%{t)\<5 i j + 



a {l + dist{i,j))v 



3.4 Weyl operators 

For operators as above we find 

W s (t) =e 1 ^^(w*'W-"A(t)) j ^ = (a; 1) ...,a;| E |,pi ) ...,p| H |) G R 2 



(75) 
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Employing the Baker-Hausdorff identity we then have, see Eq. d56*l l. 



(76) 



i.e., 



W t (t),W e 



< 



< 



UWU'W E (|c^(*)| + |c^(*)| + |c^(*)| + |c^W|) 



hi«'ii E (ii[^(*).^]ii + ii^(*).Pi]ii 

teases' 

<«ll E f(dist(ij)), 



(77) 



where 



E f(dist(i,j))= E E <W,,V- (78) 

iea,jes' d=dist(s,s') ies.jea' 

We now proceed by showing how to restrict the latter sum to subsets of 3 and S'. As 
one has to cross the boundary of a set to find a path to a site outside that set, there exist 
for all i e S, j e S' sites fc 6 <9S, Z <E &E' such that 



dist(i,j) = dist(i, k) + dist(k, I) + dist(l,j). 



(79) 



Thend = dist(i, j) requires dist(i, k) and dist(l,j) to be smaller than d— disi(H,S') =: 
r as dist(k, V) > efot(£, H'). We may thus write 



$dist(i,j),d — fidist(i,j),d> 

iea.jea' ieds r ,jeds' r 



where we denoted by 



<9A r = [J {j G A I dist(i,j) < r} 

iedA 



(80) 



(81) 



the set of lattice sites that are within A and within a layer of thickness r around the 
surface of A, for which we have 

\dA r \ < \dA\ sup \ {j £ A\dist(i,j) < r}\ 

i&dA 



< \dA\ sup \{j £L\ dist(i, j) < r} | 



(82) 



|&4| sup < |M| 1 + 



i£L 



1=0 



1 = 1 
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Hence 



(83) 



X] s dist(i,j),d < min{|<9S r | , \d3' r \} sup S d (j) 

< CD miniias.ijss;!}^- 1 

/ r ^ 

< c D mm{\dZ\ , \dE'\} d D ^ 1 + c D J] P" 1 ) 

< c D min{|d~| , |aH'|} d 15 - 1 (1 + c D r D ) . 

To summarize, 

J2 f( dist d,j)) <C2,min{|aS| > |flS'|} 

oc 

x £ /(^^(l + coCd- dis^S')) 73 ). 

(84) 

Under the assumptions of Theorem|2] e.g., we may choose 

(||p|| 11X11 \ e di - j lo s( CT M,j) 
"" + ""+ 2 -—7 iY' (85) 
\/W px W VW xp W J (l - (er/dy) 2 ) 

i.e., forer < dist(E.,E/)/R =: <fe,~', 

CD ||e||ri|min{|aS|,|55'|}( 7 ^+ J§L + 2 ) 

< , W " " x V " " V^'M-A^) 

J yfe(l-(cr/d s , s 2 ) 



x £ e dlog(eril/( (2 +di S t(H,H')))/ft (d + ^ + 

(86) 

where we have for the sum the following upper bound (z := er/ds s') 

oo 

dist D - l {~, H') £ e rfl°g(crH/(d+* St (H,H')))/« (d + ( X + Cc<i D) 



d=0 



(87) 



3=0 



4 Summary 

In this work, we have presented Lieb-Robinson bounds for harmonic lattice systems 
on general lattices, complementing and generalizing work in Refs. |0] (see also Ref. 
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1 331 13411 ). We found a stronger than exponential decay in case of local interactions, 
and an inheritance of the decay behavior in case of algebraically decaying interactions. 
For the case of the Klein-Gordon field, we found the exact locality emerging from 
the approximate locality in the Lieb-Robinson sense. Specific attention was devoted 
to the time evolution of Weyl operators, which are an important class of operators in 
harmonic lattices. As such, this work provides a framework to study non-equilibrium 
dynamics in harmonic lattice systems in a general setting. 
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